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Starting from the symmetry aspects and tight-binding method in combination with first-principles 
calculation, we systematically derive the low-energy effective Hamiltonian involving spin-orbit cou- 
pling (SOC) for silicene, which is very general because this Hamiltonian applies to not only the 
silicene itself but also the low-buckled counterparts of graphene for other group IVA elements Ge 
and Sn, as well as graphene when the structure returns to the planar geometry. The effective 
Hamitonian is the analogue to the first graphene quantum spin Hall effect (QSHE) Hamiltonian. 
Similar to graphene model, the effective SOC in low-buckled geometry opens a gap at Dirac points 
and establishes QSHE. The effective SOC actually contains first order in the atomic intrinsic SOC 
strength £o, while such leading order contribution of SOC vanishes in planar structure. Therefore, 
silicene as well as low-buckled counterparts of graphene for other group IVA elements Ge and Sn 
has much larger gap opened by effective SOC at Dirac points than graphene due to low-buckled 
geometry and larger atomic intrinsic SOC strength. Further, the more buckled is the structure, the 
greater is the gap. Therefore, QSHE can be observed in low-buckled Si, Ge, and Sn systems in an 
experimentally accessible temperature regime. In addition, the Rashba SOC in silicene is intrinsic 
due to its own low-buckled geometry, which vanishes at Dirac point K, while has nonzero value with 
k deviation from the K point. Therefore, the QSHE in silicene is robust against to the intrinsic 
Rashba SOC. 

PACS numbers: 73.43.-f, 85.75.-d, 73.22.-f, 71.70.Ej 



I. INTRODUCTION 

Silicene, as the counterpart of graphene for silicon, 
with slightly buckled honeycomb geometry has been syn- 
thesized through epitaxial growth 1 . This novel two- 
dimensional material has attracted considerable atten- 
tion both theoretically and experimentally recently, due 
to exotic electronic structure and promising applications 
in nanoelectronics as well as compatibility with current 
silicon-based electronic technology^^.The structure of 
silicene is shown in Fig. [l] In the absence of spin-orbit 
coupling(SOC), the band structure of silicene shows lin- 
ear energy spectrum crossing at the Fermi level around 
the Dirac points K and K* of the Hexagonal Brillouin 
zon e 1 2 * 3 * 5 * 6 -*, which is similar to graphene case. 

Quantum spin Hall effect (QSHE), a new quantum 
state of matter with nontrivial topological property, has 
garnered great interest in the fields of condensed mat- 
ter physics and materials science due to its scientific 
importance as a novel quantum state and the techno- 
logical applications in spintronics^Sl. This novel elec- 
tronic state with time reversal invariance is gapped in 
the bulk and conducts charge and spin in gapless edge 
states without dissipation at the sample boundaries. The 
existence of QSHE was first proposed by Kane and Mele 
in graphene in which SOC opens a band gap at Dirac 
points 10 . Subsequent works, however, showed that the 
SOC is rather weak, which is in fact the second order 
process of the atomic intrinsic spin orbit interaction for 
graphene, and the QSHE in grap hene can occur only 
at unrealistically low temperatur e! 11 * 12 *. So far, there 
is only one system, two-dimensional HgTe-CdTe quan- 
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FIG. 1. (color online). The lattice geometry of low-buckled 
silicene. (a), (b) The lattice geometry from the side view and 
top view, respectively. Note that A sublattice (red or gray) 
and B sublattice (yellow or light gray) are not coplanar. (c) 
The angle 6 is defined as being between the Si-Si bond and 
the z direction normal to the plane, (d) The relativistic band 
structure of low-buckled silicene. Inset: zooming in the energy 
dispersion near the K point and the gap induced by SOC. 

turn wells, where QSHE is demo nstra ted^ 3 ^*, in spite of 
some other theoretic suggestion s 15 * 16 !. Recently, there is 
evidence for helical edge modes in inverted InAs-GaSb 
quantum wella^ experimentally^. Nevertheless, HgTe 
quantum wells and other systems more or less have se- 
rious limitations such as toxicity, difficulty in processing 
and incompatibility with current silicon-based electronic 
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technology. Therefore, the true realization of QSHE in 
silicene is very worth while to expect. Silicene and two- 
dimensional low-buckled honeycomb structures of germa- 
nium and tin with QSHE are promising candidates for 
constructing novel spintronic devices. 

Using the first-principles method, we have recently 
demonstrated that silicene and two-dimensional low- 
buckled honeycomb structure of germanium can realize 
the QSHE by exploiting adiabatic continuity and the di- 
rect calculation of the Z 2 topological invariant!^ with 
a sizable gap opened at the Dirac points due to SOC 
and the low-buckled structures in our recent Letter^. Al- 
though the electronic structure, especial linear energy 
spectrum of silicene, at low energy is similar to that 
of graphen d 20 * 21 !, the low-buckled geometry makes the 
derivation of low-energy effective model Hamiltonian not 
as clear as graphene. Motivated by the fundamental in- 
terest associated with QSHE and SOC in silicene, we at- 
tempt to give a low-energy effective model Hamiltonian 
to capture the main physics. 

The paper is organized as follows. In Sec. II we 
briefly describe SOC in silicene from symmetry argu- 
ments. Thus, we introduce a next nearest neighbor tight- 
binding lattice model Hamiltonian to include time rever- 
sal invariant spin-orbit interaction. Section III presents 
the derivation of our low-energy effective model Hamilto- 
nian step-by-step. We investigate in detail the effective 
spin-orbit interaction including intrinsic Rashba SOC. In 
Sec. IV, a comparison of gap opened by SOC obtained 
from between our previous first-principles results and the 
current tight-binding method is made. As an application 
of our model Hamiltonian , we also study the counter- 
parts of graphene for other group IVA elements Ge and 
Sn, which are low-buckled structure according to first- 
principles calculations. We conclude in Sec. V with a 
brief discussion and summary. 

II. LATTICE MODEL HAMILTONIAN 
INCLUDING SPIN-ORBIT COUPLING IN 
SILICENE FROM SYMMETRY ASPECTS 

In general, SOC in Pauli equation can be written as, 

where V (F) is potential energy (force) , p is momentum, 
h is Plank's constant, mo is the mass of a free electron, 
c is velocity of light, a is the vector of Pauli matrices. 

For graphene as shown in Fig.[2ja), the nearest neigh- 
bor SOC is zero due to its mirror symmetry with respec- 
tive to an arbitrary bond, while the next nearest neighbor 
SOC is nonzero. According to symmetry 

H so = i"/2 (F\\ x dijj ■ a = it 2 Vijcr z , (2) 
where = d j xd j } j 2 anc ] t 2 are undetermined parame- 

| di x dj | 

ters, di and dj are two nearest bonds connecting the next 




FIG. 2. (color online). The atomic intrinsic spin-orbit in- 
teraction from symmetry aspects, (a) The nearest neighbor 
force F\ vanishes, while the next nearest neighbor force Fj| 
is non-zero in horizontal plane, (b) The next nearest neigh- 
bor non-zero force Ff equals to negative Ff in perpendicular 
direction. 

nearest neighbor dij. 

For silicene, the nearest neighbour SOC is zero, while 
the next nearest neighbor SOC is nonzero which can be 
divided into two parts, namely parallel with and per- 
pendicular to the plane, respectively, according to two 
components of electric field force (see Fig. Ill), of which 
the perpendicular component is due to A sublattice and 
B sublattice being noncoplanar. 

As the first part, the force parallel with the plane is 
taken into account. This case is similar to the graphene. 

H sol = £72 \ F\\ x dijj ■ a = it 2 Vija z . (3) 

For the second part, the force perpendicular to the plane 
is taken into account as shown in Fig. [2jb). 

H so2 = i 7l (a x 4) ' F±S„ = itiHij (a X d?) , (4) 

where 4 = dy/|dy|, 7i and t\ are undetermined param- 
eters, = ±1 for A site (B site). 

Finally, we introduce a second nearest neighbor tight- 
binding model 

(ij)a {{ij))afl 

/ — ( 5 ) 

- ih 2J (Mjc] a x 4) c jP . 

The first term is the usual nearest neighbor hopping term. 
The second and third terms are effective SOC and intrin- 
sic Rahsba SOC. The three parameters t,t 2 ,ti are given 
explicit expression forms in the following derivation by 
using tight-binding method. 

By performing Fourier transformations, we obtain the 
low-energy effective Hamitoniam around Dirac point K 
in the basis {\A), \B)} ® {t,|} 

Tjeff _ ( h n v F (k x +ik y )\ . . 

~ \v F (k x -ik v ) -h n ) ' 



3 



v F 



V3 



at, h\\ 



-3v3t2<r» 



t\Cl(kyO x k X (Jy) 



Around Dirac point K* in the basis {\A), \B)} ® {til}, 
we have 



H 



v F {k x - iky) 



v F (k x +ik y ) hn 



(7) 



The two effective Hamiltonian Eqs. ^ and ^ should be 
related by the time-reversal operation. 

From the symmetry aspects analysis, we obtain 
the effective Hamiltonian for silicene that shown as 
Eqs. ([5]) ([6]) (0. However, magnitude of the parame- 
ter in effective model and microscopic mechanism such 
as geometry enhanced effective SOC 5 ' etc are quite un- 
clear. In order to study such effect, we need to construct 
the effective Hamiltonian from the atomic tight binding 
Hamiltonian. 



III. LOW-ENERGY EFFECTIVE 
HAMILTONIAN FROM TIGHT-BINDING 
THEORY 

A. Low-energy effective Hamiltonian without SOC 

The outer shell orbitals of silicon, namely 
3s,3p x ,3p y ,3p z , are naturally taken into account in 
our analytic calculation. As shown in Fig. [I] there are 
A and B two distinct sites in the honeycomb lattice 
unit cell of silicene. Therefore, in the representation 
{\pi),\p^)M)Api)As A ),\p^)A^),\s B )} (For sim- 
plicity, the Dirac ket is then omitted over the following 
context) and at the K point, the total Hamiltonian in 
Slater Koster frame reads 



( 














V 3 ' 


-iv; 











V3 


iv; 

















V3 











-y'1 


-iV[ 


v 2 ' 
















-iv; 


v[ 


-i\ 














A 


-v% 


iV 2 ' 





V3 





-v[ 




















iVi 


v[ 


-iV 2 











V 





V-2 


w 2 











A 



\ 



We rewrite the total Hamiltonian in new basis 
{p A , s A , ^ ,pf , s B , ip A , (pa, tp 4 } 

H Q ^H 1 = UtHoUi, 



where V 1 ,V 2 ,V 3 are related to bond parameters (V ssa 
etc), the detailed derivations are shown in Appendix A. 
To diagonalize the total Hamiltonian, we take two steps. 
Firstly, we perform unitary transformation 
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H, = 



(10) 

where Ui is the unitary matrix that connects the new 
basis and original basis, Vi — 2V 1 , V 2 = \/2V 2 , V 3 — 
s/2V 3 . 

Secondly, the new Hamiltonian Hi can be separated 
to three decoupled diagonal blocks, which are named as 
Ha, Hb and He, respectively. Ha reads in the basis 
{p A ,s A ,Lp B } 



-iVa 
II t = I A iV 2 

iV 3 -iV 2 



(11) 



Its eigenvalues £1, e 2 , e 3 satisfy the eigen-equation 

E 3 — AE 2 — (V 2 + V 3 )E + AV 3 = 0. (12) 

Since the above equation is cubic equation, the eigenval- 
ues and eigenvectors of Ha can be analytically obtained. 
We perform unitary transformation {(pi, 02, ^3} = 
{p A ,s A ,(p^ } U A , where 



U A = 
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te 3 
Q3 V3 



(8) with the normalization factors 



(13) 



(A - Si)Va 



For simplify, Ua is expressed as Ua = where Uij 

is the matrix element of Ua- We rewrite Ha in the new 
basis j>i, 02, 03 } 



H A ^H' A = U\H A U A = I e 2 



£1 
£ 2 
£ 3 



(14) 



The above technique in Ha can also apply to the second 
diagonal block Hb which reads in the basis {p B , s B , <p A } . 
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Hb satisfy the same eigen-equation Eq. (12 1. Its eigen- 
values £4,65,66 satisfy: 



£4 — El, £5 — £2, £6 — ^3- 



(15) 



The eigenvectors Ub of are a little different from that 
oiH A . 



Mil Ul2 "13 
{7g = I -U21 -""22 -U23 
M31 M32 M33 



(16) 



where u%j is the matrix element in the unitary matrix Ua 



as present in Eq. ( [13] ). We define the unitary transfor- 
mation {0 4 ,0 5 ,0 6 } = {pf , s s , ^} {7 B . Obviously ff B 
is diagonal in new basis. He itself is diagonal. We define 

we have found a uni- 



97 = <^3,<P8 

From Eq. (p) to Eq. (fl6| 



tary transformation C/2 that connects the original ba- 



sis {pf,s A 



<p B ,p B , s B , ip A , (fs, ^4} and the new basis 
{0i, 04, 02, 05, 03, 06, 07, 0s}- Under such unitary trans- 
formation, Hi will be diagonal. 

Combining the above two steps, we finally find the 
new basis {0i, 04, 02, 05, 03, 06, 07, 0s} an d the unitary 
transformation matrix U — U\XJi which diagonalize the 
original Hamiltonian Hq. The results were summarized 



{01, 04, 02, 05, 03, 06, 07, 08} 



(17) 
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H = U^H a U, 
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£1 
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0000 £ 3 00 

00000 £ 3 

Vi 

\0000000 -Vi/ 

So far, the diagonal Hamiltonian has been obtained. No- 
tice that the interesting structure is low-buckled, which 
means V3 is small due to the angle 9 approaching to 90° . 
When V3 is small, the three roots of the eigen-equation 



Eq. (12) reads 
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(19) 
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Next, we determine the Fermi energy of silicene. Due 
to its half filling, there are four eigenvalues below the 
Fermi energy. According to Eqs. ( 18 ) and ( 19 1 , the eigen- 



values £3, V\ are below e\ while the others are above £1, 



so the Fermi energy locates around s\. Thus, 0i and 04 
are low-energy states which have explicit forms 



ui\p z 



u 2 is A + u 31 [^=(p B -ip*)], 
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(20) 



In order to study the low-energy physics near the Dirac 
K point, we perform the small k expansion around K by 
k — > k + K and project the Hamiltonian to the represen- 
tation {0i, ^4}. We keep the first order term of k 



H K = E 1 l 



U2 



vpk^ 
vpk- 



(21) 



with the Fermi velocity vf 

{V pjm sin 2 9 + V ppu cos 2 9) - 
2wnU2i cos9V spa - l\u 31 \ 2 sin 2 9 (V ppa 

k-^- — krp 4~ iky , k — — k'£ iky , 



u 21^ssc 
' Vppir)], 



(22) 



where a is the lattice constant, 9 is the angle between 
the Si-Si bond and z direction. Notice we have let h — 
1. So when we calculate Fermi velocity vf, should be 
considered. 



Eq. (21 I and Eq. (22) are the final low-energy effective 



Hamiltonian without SOC. The two important results 
can obtained from these two equations. Firstly, similar to 
graphene, the low-buckled silicene remains gapless with 
linear dispersion. Secondly, vp here is original from all 
parameters V P rrK,Vpp a ,V ssai V spa , while the Fermi veloc- 
ity vp is only determined by parameter V pp7r in graphene 

a ) 

2 V PP7T u ) ■ 
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B. Low-energy effective Hamiltonian with SOC 



The form of SOC Hamiltonian H so is given in the rep- 
resentation {p A ,pf,p A ,p A ,s A ,p B ,p B : s B }®{t,l} (Ap- 
pendix B). We know that Hamiltonian without SOC in 
the basis set {0i, 4 , 2 , 05, 03, 06, 07, 0s}® {f, 1} is diag- 
onal from the above depiction. The two representations 
are related by unitary transformation (Eq. (|17[)) 



U so = U®h, 



(23) 



where 7 2 is 2 x 2 identity matrix for the spin 
degree of freedom. In the representation of 
{01,04,02,05,03,06,07,0s} ® {t4}, SOC Hamil- 
tonian H~ n and total Hamiltonian H read 



H 



^so — Ug H S0 U S0 , 



H =H n ®h 



H„ 



(24) 



The first 4x4 diagonal block in SOC Hamiltonian H so 
is no other than the first order SOC, which reads at the 
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Dirac point K in the basis |( 
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(25) 



(26) 



where U31 is the corresponding matrix element in U a- 
In the following, we explain the microscopic mechanism 
leading to the above equation. The intrinsic effective first 
order SOC can be summarized as: 



(27) 



where V means the nearest neighbor direct hopping, 
£0 represents the atomic intrinsic spin-orbit interaction 
strength. The whole process can be divided to three 
steps. Take pf for example. Firstly, due to the low- 
buckled structure , pf couples to p^(see in Eq. (Ill and 
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Eq. ( 20 1 . Carrier in pf orbit directly hops to the nearest 
neighbor p B orbit. Secondly, when the atomic intrinsic 
SOC is introduced, the energy of p B will split with spin 
up carrier shifting — ^ while spin down carrier shifting 
%. The third, carrier in p B directly hops to another 
nearest neighbor pf orbit. The SOC process in the pf is 
analogous to that of pf except that pf orbit couples to 
the pf orbit. The difference leads to the opposite mag- 
nitude of effective SOC . During the whole SOC process, 
the atomic intrinsic SOC takes effect for only one time. 
Therefore, the effective SOC is proportional to Co- A 
brief sketch of the process is shown in Fig. |3jh). 

For most of time, we focus on the low-buckled geom- 
etry with small V 3 . According to the Eqs. (p| (fl9]| (po) , 
Al** reads 
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Especially, when low-buckled geometry returns to planar 
structure such as graphene^ = 90°), the above formula 
becomes A^* = and the first order SOC vanishes. Phys- 
ically, when 9 = 90°, pf orbit is orthogonal with p% and 
Py orbits. Therefore, the directly hoping from pf to p B 
is completely forbidden. The SOC process described in 



Eq. (27) cannot happen. 



The effective second order spin-orbit interaction is also 
deduced, whose detail derivation is described in Ap- 
pendix C. Here, we do not intend to repeat them again 
but just to quote some expressions there Eq. (|C9| 



HIT = ~ H n (H a - ex) Hi, 



(29) 



where H n take from the first row to the fourth row and 
the fifth column to the sixteenth column of the above H so , 
H a is the direct product matrix between the lower right 
6x6 diagonal matrix of ff (Eq. ( |l8| ) and 2x2 identity 
matrix, E\ is the eigenvalue aforementioned. When V3 is 
small, the effective second order SOC Hamiltonian reads 
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We analyze the microscopic mechanism for the \ 2 s f . 
Without SOC, the low-energy H n and the high energy 
H a are decoupled. However, in the presence of the atomic 
intrinsic SOC, H n and H a are coupled together. A detail 
analysis is shown that A 2 *' can be summarized as the 
process as: 
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where V means the nearest neighbor direct hopping, 
Co represents the atomic intrinsic spin-orbit interaction 
strength. During the process, the atomic intrinsic SOC 
takes effect for twice. Thus, the effective SOC is the sec- 
ond order of Co- A brief sketch of the process is shown in 
Fig. |3jc). We notice that in graphene, the second order 
A 2 ** is the leading order of eff ective SOC and had been 
studied in reference d 11 | 12 | 21 | 2 2J . 



C. Intrinsic Rashba SOC in silicene 

The extrinsic Rashba SOC in graphene is due to a per- 
pendicular electronic field or interaction with a substrate 
which breaks the mirror symmetry, while the intrinsic 
Rashba SOC in silicene is due to its own low-buckled ge- 
ometry. Around the K point, the Hamiltonian contain- 
ing k deviation from the K point in the representation 
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where Hq is given in Eq. 
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Through the unitary transformation matrix U (Eq. ([T7|)), 
in the representation 04, <j> 2 , <f>5, $3, 4>q, 4>r, (f>s}, we 
have 

Hq (fc) — ► ^ (A) = E^£f (*) ^ = ^0 (*) + #o, (35) 



where _ff is given in Eq. (18 1. We mainly focus on the 



terms containing k deviation from the K point 

H (k) — > H (k) = H'q (k) <E>I 2 + H so , (36) 



where H so is given in Eq. (24 1. According to the Eq. ( |C9| , 
the total second order Hamiltonian reads 

H' eff (k) = -H non (k) (H a - e^- 1 Hl n (k) , (37) 

where H non (k) takes from the first row to the fourth row 
and the fifth column to the sixteenth column of H (fc). 
The Hamiltonian H e j^{k) can be divide into two parts 



H eff (k)=H^ + H R (k) 



(38) 



where if 2 ^' is given in Eq. (30). HrQc) is intrinsic 
Rashba SOC in silicone, which can be written as around 
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where the purely real Xr reads 
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From the above equations, we know that Hr (fc) is ex- 
actly zero at Dirac point K , while i/^ (fc) has nonzero 
value with k deviation from the K point. Moreover, 
when the structure returns to the planar structure, 9 — 
90°, Aii = 0, the intrinsic Rashba SOC vanishes even 
when k deviating from K. Therefore the intrinsic Rashba 
is entirely caused by the low-buckled geometry. The in- 
trinsic Rashba SOC is quite different from the extrinsic 
Rashba SOC, which arising from a perpendicular elec- 
tronic field or interaction with a substrate leading to mir- 
ror symmetry broken in some direction, has finite mag- 
nitude at Dirac point K. 



IV. RESULT AND DISCUSSION 



Finally, in combination with Eqs. (|21J) (|25| ([30j (|39| , 
we obtain the entire low-energy effective Hamiltonian 
around Dirac K acting on the low-energy states </>i and 
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where I 2 is 2 x 2 identity matrix for the spin degree of 
freedom, I4 is 4 x 4 identity matrix and A so = A™* + 
A 2 q', vf is given in Eq. (22). Through the time-reversal 



operation, the entire Low-energy effective Hamiltonian 
around Dirac K* reads 



H e K { f (e) = (e 1 -XH t )h + 



( —/in Vpk- 
y v F k+ hu 



(42) 



The effective Hamiltonian deduced from atomic tight- 
binding method has the similar formulas as that from 
the symmetry aspects. Comparing Eq. (411 and Eq. 
we obtain 



2\/3' 



3 a 



A S o 



(43) 



The above parameters t,t2,ti are undetermined in the 
second nearest neighbor tight-binding model Eqs. ([5]) and 
^ from the above symmetry analysis. Here, in combi- 
nation with Eqs. (l22f p8l (pll pOl) (El , we can not only 



give their explicit expressions, but also specify the mag- 
nitudes of the three parameters through vp, A so (A so = 
Ago* + As* 4 ), and \r, whose values are presented in Table 

in 

In the following, we discuss the physic meanings of 
our obtained low-energy effective Hamiltonian. First of 
all, the low-energy effective Hamiltonian is analogous 
to the first QSHE proposal in graphene except the in- 
trinsic Rashba SOC term H R (fc) 19 . The SOC inducing 
mass term to Hamiltonian opens gap at the Dirac points. 
Moreover, from K to K* the mass term changes its sign 
and the band is inverted. Therefore, the low-buckled 
silicene is also the QSHE system. The QSHE can be 
observed experimentally when the Fermi energy locates 
inside the gap and the temperature is below the minimal 
energy gap. The existence of the QSHE in silicene has 
been studied in our recent work using the first-principles 
method combined with the direct Zi calculation^!. 

Secondly, the energy gap in low-buckled silicene is 



much larger than that in graphene. The Eq. (41 1 results 
in a spectrum E(k) = ±\/ (vp + a 2 X 2 R ) k 2 + \ 2 so . There- 
fore, the energy gap is 2A so at the Dirac points. Due to 
the low-buckled geometry, not only the second order SOC 
A^o* but also the much larger first order SOC Aj** exist. 
In Fig. [3] we show the variation of gap with the angle 0. 
With 9 deviating from 90°, the gap induced by A^' for 
silicene is nearly unchanged while the gap induced by A**' 
increases rapidly. The larger is the angle, the greater is 
the gap. Especially, the gap can reach to several meV for 
just little buckled, therefore the QSHE can be observed 
in an experimental observable temperature regime. 

Thirdly, due to the low-buckled geometry, the effec- 
tive Hamiltonian also contains the intrinsic Rashba SOC 
term. Such term leads to interesting properties. On the 
one hand, since it vanishes at Dirac point, the minimal 
bulk energy gap 2A so will not be affected by the intrinsic 
Rashba SOC. Therefore, it does not diminish the temper- 
ature window for experimentally observing the QSHE in 
silicene. On the other hand, due to the nonzero values 
of the Rashba SOC term, spin is not a good quantum 
number. Thus, the spin Hall conductance is no longer 
quantized in silicene. The intrinsic Rashba SOC is en- 
tirely different from the extrinsic Rashba SOC, which has 
finite value at Dirac point K, can destroy the QSHE 19 . 
The presence of intrinsic Rashba SOC may provide a way 
to manipulate the spin in silicene without destroying its 
QSHE state. 



> 





90 95 100 105 110 

6 (degree) 

FIG. 3. (color online), (a) The variation of gap opened by 
SOC at Dirac point with the angle 9 for silicene. Black line 
marks the total gap. Red or gray line and dot line mean gaps 
opened by the first order and second order SOC respectively, 
(b), (c) Brief sketches of the intrinsic effective first and second 
order spin orbit interaction. 



Fourthly, the entire low-energy effective Hamiltonian 
applies to not only the silicene itself but also the low- 
buckled counterparts of graphene for other group IVA 
elements Ge and Sn, as well as graphene with the planar 
geometry. These different structures correspond to the 
different angles 8. Therefore, in this sense, the effective 
Hamiltonian is quite general. 

In Table [i] A 1 *', A^*, Xr, the gap, and vp caused 
by SOC at Dirac point K in graphene, silicene, gelicene 
and snlicene, similarly corresponding to two-dimensional 
low-buckled Ge and Sn, are obtained from tight-binding 
meth od b y using the typical parameters values from Ta- 
Notice that A^* is slightly larger than A 



III 



Is/ 



ble 

due to the huge SOC strength with the magnitude of 
eV in snlicene, while A^* is much smaller than A 1 *' in 
other systems. For comparison, from the first-principles 
method, we present the corresponding gaps too, which 
agree with our tight-binding method results in order of 
magnitude. We also give the carrier Fermi velocity Vf 
around Dirac point K from first-principles and the cur- 
rent tight-binding method. Since we only focus on the 
low-buckled geometry, our calculation shows that the car- 
rier Fermi velocity does not significantly change with 6. 

Notice that those bond parameters presented in Ta- 
ble |III| used in Table [I] come from the corresponding di- 
amond structure(sp3 hybridization) actually except the 
graphene(sp2) case. However, considering that low- 
bulked structures are more closely to sp2 hybridization 
and bond parameters of sp2 hybridization will be a little 
different from those of sp3 hybridization, thus, through 
slight improvement of these bond parameters we expect 
that the tight-binding gap can better match with the 
first-principles results. 
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TABLE I. An application of the low-energy effective Hamiltonian (Eq. (41 1). The terms of this general low-energy effective 
Hamiltonian for these different systems corresponding to the different angles 9 of the lowest energy geometry are given. The 
lattice constants a with unit A and angles 8 for the lowest energy geometry are obtained from the first-principles calculation. 



A 



A: 



and Xn caused by low-buckled geometry and SOC, with units meV at the Dirac point K are presented obtained 



from our tight-binding model using hopping parameters in Table [nil The gap with unit meV opened by SOC at Dirac point 
K from first-principles and the current tight-binding method is presented. We also give the carrier Fermi velocity Vf with unit 
10 5 m/s around Dirac point K from first-principles and the current tight-binding method. 



System 


a(A) 


e 


Xlo (meV) 


\2st 


X R 


Gap (meV) (TB) 


Gap(FP) 


v F (I0 5 m/s) (TB) 


vf(FP) 


graphene 


2.46 


90° 





1.3 xl(T 3 





2.6xl0" 3 


0.8x10" 




9.80 


8.46 


silicene 


3.86 


101.7° 


3.9 


7.3xl0~ 2 


0.7 


7.9 


1.5E 


b 




5.52 


5.42 


ge(licene) 


4.02 


106.5° 


43 


3.3 


10.7 


93 


23. £ 


b 




4.57 


5.24 


sn(licene) 


4.70 


107.1° 


29.9 


34.5 


9.5 


129 


73.5 




4.85 


4.70 



a Reference II II 
b Reference \E\ 



V. SUMMARY 

In summary, based on the symmetry aspects and tight- 
binding method combined with first-principles calcula- 
tion, we derived the low-energy effective Hamiltonian for 
silicene, which is very general because this Hamiltonian 
applies to not only the silicene itself but also the low- 
buckled counterparts of graphene for other group IVA 
element Ge and Sn, as well as graphene when the struc- 
ture returns to the planar geometry. The low-energy ef- 
fective Hamiltonian is indeed QSHE with its form similar 
to Kane-Mele's first graphene QSHE Hamiltonian except 
the intrinsic Rashba SOC term H_n(k). However, the 
effective SOC in low-buckled geometry is actually first 
order in the atomic intrinsic SOC strength £o> while the 
planar structure in graphene leads to the vanishing of 
the leading-order contribution. Therefore, silicene as well 
as low-buckled counterparts of graphene for other group 
IVA elements Ge and Sn has much larger gap opened 
by effective SOC at Dirac point than graphene due to 
low-buckled geometry and larger atomic intrinsic SOC 
strength. Further, the larger is the angle, the greater 
is the gap. Therefore, QSHE can be observed in an ex- 
perimentally accessible low temperature regime in these 
low-buckled systems. In addition, the Rashba SOC in 
silicene is intrinsic due to its own low-buckled geometry, 
which vanishes at Dirac point K, while has nonzero value 
with k deviation from the K point. As a result, though 
spin Hall conductance is not quantized, the QSHE in the 
silicene is robust against to such intrinsic Rashba SOC. 
This is entirely different from the extrinsic Rashba SOC 
due to a perpendicular electronic field or interaction with 
a substrate, which is independent of k, has finite value 
at Dirac points, and can destroy the QSHE. 
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Appendix A: Ho Matrix 

In the representation {pf,pf,Py,p^,s A ,Py,p^,s B ] 
the total Hamiltonian reads 



Ho 





(Al) 



Here, H T and are 2x2 and 6x6 matrices, respec- 
tively. The non-diagonal block H n coupling and H CT 
is 2 x 6 matrix. In the following derivation, the energy 
level of 3p orbital is set as energy zero point. The matrix 
E describes the on-site energy of different atomic orbitals, 
which can be written as 



(A2) 



where A is the energy difference between the 3s and 3p 
orbitals. Actually, here we have assumed these basis are 
orthogonal when centered on different sites for the sake 
of simplicity. We choose the coordinate system in which 
the unit cell has primitive vectors 





(o 








E = 


















A 
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The lattice constant a is defined as the nearest distance 
of lattice point at the same sublattice, which is 3.86 A for 
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silicene from our first principles calculation^. The three 
nearest neighbor translation vectors are 



r a ( V3 1 

dl = 

r a f y/3 1 

d-2 = , - , cot I 

^ 2 '2' 

d 3 = (O,-l,cot0) . 
v3 



(A4) 



As shown in Fig. [I] the angle 9 is defined as being between 
the Si-Si bond and the z direction normal to the plane. 
The corresponding reciprocal lattice vectors are 




V3 



(A5) 



The Dirac point K is chosen to be K = | (b\ — b%j = 

(|^,0), as well as K*=— K. The matrix T describes the 
hopping between two sublattices, which is given in Ta- 
ble |n] by the Slater-Koster formulfP. In Table |nj the 
four bond parameters V ssa , V sptT , V ppcr and V VV7I corre- 
spond to the cr and 7r bond formed by 3s and 3p or- 
bitals, whose numerical values given in Table |ITT| specify 
our model quantitatively. The hopping matrix elements 
in the momentum space read 



TABLE II. The matrix elements for the nearest neighbor hop- 
ping between s and p orbitals are considered as functions of 
the direction cosine I, m and n of the vector from the left 
orbital to the right orbital. Other matrix elements are found 
by permuting indices. 



is,S VsscT tx,X I Vpp(j^(l-^ ) Vpp-jr 

ts,X IVspa tx^y ltTl( t Vppa--Vpp7T^) 

ix : s ~lVspa ty,z TTLTl ( Vppcr " Vpprr ) 



TABLE III. The magnitudes of hopping parameters. The 
energy units are eV. The strength of SOC £o is obtained 
from first-principles calculation except Sn. 



System 
Graphene 

Silicene 
Ge(licene) 
Sn(licene) 



Vss(T 

-6.769 
-1.93 
-1.79 



5.580 
2.54 
2.36 



5.037 
4.47 
4.15 



Vpp7T 

-3.033 
-1.12 
-1.04 



A fr> 
-8.86^ 9x10" 

-7.0: 

-8.0! 



-2.6245 2.6504 1.4926 -0.7877 -6.233 




a Reference 24 
b Reference 1251 
c Reference [Til 
d Reference \E\ 
e Reference 1261 
f Reference [271 



Appendix B: H so Matrix 



t(k) =Y^t{di)t 



ik ■ di 



(A6) 



Therefore, the matrix T and H n at the Dirac point K 
can be written as 



-V x -iV x v 2 
-iV x V x 



-iV 



-Vo iVn 



H = ( F 3 ~iV 3 
n I V 3 ' 1V3 



V 



v 



- Sin 2 9 (Vppir — Vppcr) , 

3 

- sm6V spa , 
3 

- sin 9 cos 9 (V pp7r - V ppcr ) 



The matrix at the Dirac point K reads 



H„ 



Consequently, Hamiltonian Hq is obtained. 




(A7) 



(A8) 



(A9) 



When in center field, Eq. ([I]) reads 

H so = £,oL ■ s. 
The above equation can also be written as 
_ (L + s_+L_s. 

n so — to 



(Bl) 



(B2) 



where s± — s x ± is y denote the plus(minus) operator 
for spin and L± = L x ± iL y denote the plus(minus) op- 
erator for the angular momentum in the selected basis. 
The SOC on the same atom is taken into account. The 
concrete SOC term can be obtained by calculating the 
mean value of the Eq. (B2|. For example, the SOC term 
between \p z f) and \p x |) reads (p z f \H so \p x I) = -f 1 
etcP^. During the derivation we may take advantage of 
the following expressions 

L + \l,m) = [I (I + 1) - to (to + 1)] 1/2 \l,m + 1) , 

L-\l,m) = [Z(; + 1)-?ti(to-1)] 1/2 |2,to-1), ( B3 ) 
L z \l, to) = m\l, to), 

where /, to represent the azimuthal quantum number and 
magnetic quantum number, respectively. A straightfor- 
ward calculation leads to the on-site SOC in the repre- 



sentation {pf,p^p y i ,p^,s A ,p^,p^ ,s B } ®{t,|} 



£0 



x h s 



(B4) 
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TABLE IV. The values of SOC among atomic orbitals that 
used in h so . A,B denote the two distinct sites. The nonzero 
SOC terms only exist in the same site. a x ,y,z are Pauli ma- 
trixes acting on the spin space. O denotes the zero matrix. 





AIB 
Pz 


AIB 
Py 


AIB 
Px 


S A/B 


AIB 
Pz 





ia x 


-Iffy 





AIB 
Py 


-ia x 





y<T z 





AIB 
px 


ia-y 


-icr z 








S A/B 





O 


o 






All elements in h qn can be found in the Table \W\ 



Appendix C: the Second order effective Hamiltonian 

In general, the Hamiltonian reads 

(CI) 



Hit H n 

Hi H„ 



We focus on the case: (i) the eigenvalues of H v is around 
energy e while the eigenvalues of H a is far away from e 
(ii) the energy scale of the non-diagonal block H n is much 
smaller than the eigenvalue value difference between H^ 
and H a . The effective Hamiltonian around energy e (or 
second order effective Hamiltonian for H^) can be ob- 
tained by the following method 2 ^. H can be rewritten 
as 



= el 



(C2) 



H-rr 



1L 



f: 





For simplicity, we omit the unitary matrix / in the above 
and the following derivation. In order to obtain the ef- 
fective Hamiltonian, one may perform a canonical trans- 
formation: 

H — > H s = er s Ue s , 



S = 



M 
-Aft 



(C3) 



where the matrix M is determined by 

[H ,S]+H„ on = (C4) 

Through simple algebraic derivation, we have 

Off, - e)M - M(H a -e)+H n = 0. (C5) 

Therefore, we can find a recursive expression for M 

M = [H n + (Hr - e)M](H a - e)- 1 
= H n {H a - e)- 1 + (H n - e)H n {H a - e)" 2 + • • • . 

We know that in silicene the eigenvalues of H a — e deter- 
mined by the energy of H a separated from those of 
are of order eV near the K point, while the energy scale 
of H n — e is nearly zero as well as H n is of order meV 
for SOC. Therefore the above recursive expression can be 
written as: 



M*H n (H a -e)- 1 . 



(C7) 



The transformed Hamiltonian has the following approx- 
imate form 



Hs = e- s Ue s = \ 



[H,5] + i[[H,5],5] 



= el 



1, 



no "t- rpH non , S] 



(C8) 



Up to the second order, the final effective Hamiltonian 
for Hit can be written as 



H eSS ~H^- X - (H n M* + MHl) 
— H w — H n (Ha- — e) 1 H} 1 . 



(C9) 
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